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We study the next-to-leading order perturbative QCD corrections to the transverse momentum- 
weighted Sivers asymmetry in semi-inclusive hadron production in lepton-proton deep inelastic scat- 
tering. The corresponding differential cross section is evaluated as a convolution of a twist-three 

■ quark-gluon correlation function, often referred to as Qiu-Sterman function, the usual unpolarized 
fragmentation function, and a hard coefficient function. By studying the coUinear divergence struc- 

. ture, we identify the evolution kernel for the Qiu-Sterman function. The hard coefficient function, 

■ which is finite and free of any divergence, is evaluated at one-loop order. 

u 

In recent years, transverse spin physics has attracted tremendous attention from both the experimental and theo- 
retical communities. As transverse spin can correlate with the transverse momentum of the partons inside a polarized 
^ ■ proton, transverse spin observables, such as single transverse spin asymmetries (SSAs), are sensitive probes of par- 
i ' ton's transverse motion and a path to three-dimensional proton tomography [H. Significant theoretical progress has 
been made in studying the single transverse spin asymmetries in the past several years. Two QCD mechanisms for 
^— I ' generating SSAs have been proposed and applied extensively to phenomenology: the transverse momentum dependent 
. (TMD) factorization approach |3J6| and the collinear twist-3 factorization approach They were shown to be 

■ closely related and, thus, provide a unified picture for the SSAs p^ . 

^.^g [ One of the most studied asymmetries is the so-called Sivers effect [l^. At the partonic level, it corresponds to an 

■ azimuthal correlation ^ S± ■ {P x k±), with S± and P being the spin and momentum vector of the polarized proton 
• ' and k± being the transverse momentum of the parton. Such correlation is encoded in the so-called Sivers function, 

. if one uses the TMD factorization formalism; or the twist-3 quark-gluon correlation function, often referred to as 
' Qiu-Sterman function, within the collinear twist-3 factorization formalism. The evolution equations for either the 
. Sivers function (l3 - [T6| or the Qiu-Sterman function [13-113] have been derived recently, which enhances the accuracy 
• • ■ of the phenomenological applications. 

A natural step forward, as a follow-up to the derivation of evolution equations, will be the computation of the next- 
to-leading order (NLO) corrections to the transverse spin-dependent cross sections. Although tremendous progress has 
I been made in the evaluation of NLO perturbative QCD (pQCD) corrections to the spin-averaged cross sections, similar 
efforts on the transverse spin-dependent cross sections are still rather limited. This is largely due to the complexity of 
such type of calculations. So far, the only NLO correction in this direction is performed for Drell-Yan production [191. 
NLO corrections to other processes will provide process-dependent corrections to the hard-part coefRcient functions, 
and can also be used to extract the universal behavior of the evolution kernel for the relevant spin-dependent parton 
distributions and/or fragmentation functions. A NLO calculation for a particular physical process, thus, provides a 
direct test of QCD factorization for the associated observables. 

In this paper we follow Ref. [l^ on the Drell-Yan production and perform a NLO calculation for the P/ij^-weighted 
Sivers asymmetry in semi- inclusive deep inelastic scattering (SIDIS). Here, Ph± is the transverse momentum of the 
final-state hadron. Since the transverse momentum is being integrated out, our result is presented within the collinear 
factorization formalism in terms of twist-3 Qiu-Sterman function, NLO hard-part coefficient function, and the usual 
unpolarized fragmentation function, as we demonstrate in detail below. The rest of our paper is organized as follows. 
In Sec. II we introduce the notation for the semi-inclusive hadron production in deep inelastic scattering and present 
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the P;ij^-weighted Sivers asymmetry at leading order. In Sec. Ill we present the NLO pQCD corrections. We first 
give the result for virtual corrections and then study the real corrections. We then combine the real and virtual parts 
to obtain the final cross section. We show that all the soft divergences cancel out between real and virtual corrections. 
The remaining collinear divergence is absorbed in the redefinition of the unpolarized fragmentation function, and the 
twist-3 Qiu-Sterman function. This provides an alternative way to derive the evolution equation for the Qiu-Sterman 
function. We conclude our paper in Sec. IV. 



II. TRANSVERSE MOMENTUM- WEIGHTED SIVERS ASYMMETRY AT LEADING ORDER 



We start this section by specifying our notation and the kinematics of SIDIS. We consider the scattering of an 
unpolarized lepton e with momentum on a transversely polarized proton p with momentum P and transverse spin 
vector Sa_, 



e{l) + p(P, 5i) ^ e(£') + h{Ph) + X, 



(1) 



where h represents the observed final-state meson with momentum Ph. In the approximation of one-photon exchange, 



we define the virtual photon momentum q 
defined as follows: 



S^{P + lf^ 



xb 



I — I' and its invariant mass 

y 



2P-q " P-l xbS' 

The differential cross section that includes the so-called Sivers effect, the sm{iph 
following form [1^ 

das ivers 



Zh 



dxBdydzh<PPhA 



Co 



Fuu + sin( 



\ psin(0,,-0s) 
J^UT 



. The usual SIDIS variables are 

) module, can be written as the 

(3) 



where (Tq = ^^(^2™ ^ ^uu and F^^'^'^ '^^^ are the spin-averaged and transverse spin-dependent structure 

functions, respectively. (j)h and ips are the azimuthal angles for the final-state hadron momentum Ph and spin vector 
Si_. In this paper, we define all our angles in the so-called hadron frame (26j . The spin-averaged differential cross 
section da / dxBdydzhd^ Ph± = aoFuu, and the P^j^-integrated spin-averaged cross section is defined as 



da 



d^PhA 



da 



dxBdydzh J ' dxBdydzhd^Ph^' 
At the same time, the transverse spin-dependent differential cross section dAa{S±_) / dx Bdydzhd'^ Phj 



(4) 



(To sm( 



5)P^" ' and the transverse momentum- weighted transverse spin-dependent cross section is given by 



d{Ph±AaiS^)) 



dAa{Si_) 



(pPf^j_£°'i^ S°' P^ 

^ dxBdydzhd'^PhA 



dxBdydzfi 

where e"^ is a two-dimensional anti-symmetric tensor with e^^ — 1, and e"^S'"P,^j^ = Ph±sm{ 
Pft,_L-weighted Sivers asymmetry is given by 



A^^^-^'''^^^~'^^\xB,y,z,) 



_ d{Ph±Aa{S±)) 
dxsdydzh , 



da 



dxsdydzh 



(5) 

Thus, the 
(6) 



There are multiple studies on the higher order pQCD corrections to the spin-averaged cross sections, see, for example, 
Rcf. [27J. As a warm-up excercise, we also calculate this cross section to NLO order, and our findings are consistent 
with Ref. [iSl- Since this result is well-known, we only give the final expression here (in the MS scheme) 

da 



dx B dydzh 



E2 [ dx dz / 2 \ / 
eg / q{x,n )Dq^h{z, 
J X z 



^i')S{l - x)dil - z) + ao^ 



"-^s \ ^ 2 



X < In 



[Pqgix)S{l ~Z)+ Pggiz)S{l - X)] + Cf 



1 + 



/dx dz ^ r, 
}Dg^h(z,fi ) 
X z 

-I 

8(5(1 -i)(5(l - z) 



+5(1-z)Cf 
+5{1-x)Cf 



(1 + ^^) 



ln(l 



2, , ln(l - z) 




(7) 
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where x = xb/x and z = Zh/z and Pqq{x) is the usual spin-averaged quark-to-quark sphtting kernel 



Pqq{x)=Cf 



+ 5{l-x) 



(8) 



Let us now concentrate on the NLO correction to the transverse spin-dependent differential cross section. This will 
allow us to study the Ph±-weighted Sivers asymmetry, as in Eq. ([B]), at next-to-leading order. We will start with 
the leading order calculation for the transverse spin-dependent cross section. At this order, the final-state hadron is 
produced through the hadronization of the quark, which comes from the virtual-photon quark scattering. In order 
to obtain a non-vanishing P/ii^-weighted transverse spin-dependent cross section (P/j^ Act(5j_)), we have to include 
the final-state multiple interactions, as shown in Fig. [TJ to provide the required phase 8]. We work in the covariant 




FIG. 1. Leading order Feynman diagrams. Left: gluon to the left of the cut. Right: gluon to the right of the cut. 
gauge. For the gluon on the left of the t — oo cut, shown in Fig. (Uleft), we have the P/^j, -weighted cross section as 



dx B dydzh 



Fig. [Illoft) 



oc j d^Ph^e'^^SlPH^ j dzDq^h{z) j dxidx2d^k^Tq^A{xi,X2,k^) 

X i-g^'' H^,ixi,X2,k^)]6'' [Ph^ - zk^) . (9) 



Here (and throughout the paper), for simplicity we only consider the so-called metric contribution |27H29l |. This means 
that we contract our hadronic tensor iJ^i, with —g^'^. The twist-3 correlation function T^^^(xi, a;2, is defined as 

^ [ dy^ d'^y± ixiP+v7 t(x2-xi)p+v- zk^-y^'^ ^ 



TqAx,,X2,k^)^ j ^^^e"^^"''^e^(-^--^^^"«=e^'=-''--(P5|V,(0)7+A+(y2-,y^)^,(2;r)|P5). (10) 

One can take advantage of 5^ {Ph± — zk±) to integrate out d^Ph±, thus P^^ = zk^. We then use k^ to convert 
to the P"*"^ field strength through integration by parts ^3^. At the same time, one realizes that the Feynman diagram 
with the gluon to the right of the cut (Fig. [Ijright)) gives no contribution to the P;j^-weighted Sivers asymmetry. 
This is because for this diagram the associated 5-function becomes 5^{Ph±), i.e. Ph± = 0, and the Pft_L-weighted 
asymmetry vanishes. This conclusion also holds true in the virtual diagram calculation. 

The required phase to generate a Sivers asymmetry comes from a pole in the propagator, which is represented by 
a short-bar in Fig. [TJ 

1 11 1 / ^c. N 
o = > (—iTTjdlXi — X2)- (11) 

{pc - ix2 - xi)Pf + le 2P-pcXi-X2+te 2P ■ ' ^ 

With this phase, we have 



5. 



J d^k^ie^PSlklTq^A{xi,X2.k^) = l.Tq^p{x,,x2), (12) 



where Tq^p{xi,X2) is the well-known Qiu-Sternian function, with the following definition: 

TqA^i,X2)= J ^^^e*^^^+^^"+^(^^-^^)^^^^"i(P5|^,(0)7+e"^^j.oP+^(y2")^,(2/r)l^5)- (13) 

Finally, the P/ij^-weighted Sivers asymmetry at LO has the following form ^ 

d{Pk.±A(T{S±)) zhdosr^ 2 fdxdz 

dxBdydz, =-^Y^] --T,,FMD,^n{z)6{l-x)5{l-z), (14) 
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where we recall that x = Xb/x and z = Zh/ z. Since we will use dimensional regularization for our NLO calculation in 
the next section, we also need the LO result in n = 4 — 2e dimension. We find that the only change is the appearance 
of f — e in the overall normalization of tio, i-e. in n = 4 — 2e dimension we have in Eq. (|14|) defined as 



-^0 = ^2 (1-^)- (15) 



III. TRANSVERSE MOMENTUM- WEIGHTED SIVERS ASYMMETRY AT NLO 



In this section we present the NLO pQCD corrections to the transverse spin-dependent differential cross section. We 
first give the result for virtual corrections, and then study the real corrections. We then combine the real and virtual 
corrections to obtain the final expression. We show that all the soft divergences cancel out between real and virtual 
diagrams. The remaining collinear divergence can be absorbed by the redefinition of the unpolarized fragmentation 
function, and the twist-3 Qiu-Sterman function. This provides an alternative way to derive the evolution equation 
for the Qiu-Sterman function. 



A. Virtual corrections 




FIG. 2. The generic Feynman diagrams for the virtual corrections to the Phx-weighted cross section. 

We first study the virtual corrections. The relevant generic Feynman diagrams are shown in Fig. [2l Here, we only 
include the diagrams which have the gluon attached to the left of the cut. This is because the diagrams with the gluon 
to the right of the cut, just like in the LO calculation, give no contribution to the P/i_L-weighted Sivers asymmetry 
because of the same (5-function S'^{Ph±)- The blob in Fig.[2jleft) is given by Fig.[3l These diagrams are pretty easy to 
cornpute since they are the same as the usual virtual corrections in the unpolarized cross section. The result is given 
by M 




FIG. 3. One-loop virtual corrections for the P;ix-weighted cross section: shown here are the corrections to the quark-photon- 
quark vertex, corresponding to the blob in Fig. [SJleft). 

On the other hand, the blob in Fig.[2Jright) is much more complicated and the explicit diagrams are given in Fig. 3] 
The calculation is lengthy and contains significant amount of tensor reduction and integration. The diagrams contain 
three types of color factors: (a) and (e) have color factors Cp; (b), (c) and (f) have color factors — l/2A*'c = Cp — Nc/2; 
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(e) 



(/) 



FIG. 4. One-loop virtual corrections for the Phx-weighted cross section: shown here are the corrections to the quark-photon- 
quark vertex with gluon attachment, corresponding to the blob in Fig. [2j right). 



(d) and (g) have color factors Nc/2. We find that the terms associated with Nc/2 cancel out and only the color structure 
proportional to Cf remains. The final result is 



d{Ph±AaiS^)) 



dxBdydzh 



Fig.Hright) 2 4tt ^ 



dx dz 



1 



Tq^F{x,x)Dq^h{z)5{l - x)6{\ - z) 



Thus, the virtual correction is given by the sum of both diagrams in Fig. [2] 



d{Ph^^a{S^_)) 



dxBdydzh 



virtual 



xCf 



2 2tt 



— —Tq^F{x,x)Dq^h 



(z)5{l - x)5{l - z) 



1 



V{l-e) 



(17) 



(18) 



B. Real corrections 



Let us now study the real corrections. In this case, we need to perform the usual fc^-expansion (also referred to 
as coUinear expansion). The techniques for A:^-expansion are well established in the literature, see e.g. Refs. 0-[ll[. 
The P/j^-weighted cross section can be written as follows: 



d{Ph^^a{Si_)) 



dxBdydzh 



Real 



j d^Ph±e"^SZPl'^^J dzDg^hiz) J dxidx2d^k^Tq,A{xi,X2,k^_) 



xkP^^[-9^''H,,{x^,X2,k^)] 



(19) 



Again, we need a phase to generate the Sivers asymmetry, which also comes from the pole in the propagators. For 
the SIDIS process we can have both soft-pole and hard-pole contributions. Soft-pole contributions come from the 
Feynnian diagrams shown in Fig. [SJ with the soft-pole marked by a short bar in the diagram. For gluon to the left of 
the cut, it arises from 



1 



1 



1 



{pc - {x2 - xi)P - ki_) +ie 2P • pc xi - X2 - I'l ■ k±_ 
where the Mandelstam variables are defined as 

S = {xP + qf, i^{p,~qf, u = {xP-pcf 



{—iT:)5{xi - X2 + vi ■ k±), 



(20) 



(21) 



Here, Pc is the momentum of the final-state quark which fragments into the observed hadron and v'^ = 2xp^/u. When 
kj_ 0, Xi = X2, i.e. the attached gluon momentum becomes zero. This clarifies the name "soft-pole" contribution. 



FIG. 5. Feynman diagrams for soft-pole contributions to the Phx-weighted transverse spin-dependent cross section. The short 
bars indicate the propagators where the soft pole arises. The "mirror" diagrams (those with gluon attached to the right of the 
cut) are not shown, but are included in the calculations. 



On the other hand, the hard-pole contributions come from the Feynman diagrams shown in Fig. [6l with the 
hard-pole marked by a short-bar in the diagram. For the gluon to the left of the cut, the phase arises as follows: 

1 1 1 



{xiP + q)^+ie 2P -qxi-XB+ie s + 



{~i7T)S{xi ~ xb), (22) 



which sets xi ~ xb = xQ"^ /{s + Q^). On the other hand, the on-shell condition for the unobserved gluon leads to 
another i5-function, 5{x2 — x — V2 ■ k±) with V2 = —2xp'^/t. Thus when k± — > 0, the attached gluon momentum 
X2 — xi = X — Xb — xs/(s + Q^), which is finite. This clarifies the name "hard-pole" contribution. 




FIG. 6. Feynman diagrams for hard-pole contributions to the P^x-weighted transverse spin-dependent cross section. The short 
bars indicate the propagators where the hard pole arises. The "mirror" diagrams (those with gluon attached to the right of 
the cut) are not shown, but are included in the calculations. 



The coUinear expansion for both soft-pole and hard-pole contributions are straightforward. After such an expansion, 
we will have a linear -dependence: 



d 



[-g'^''Hf_,^{xi,X2,k±)] ocp^ 



pp 



This linear P^_|_ will combine with P^^ in Eq. (|19l) to become P^j^, 

p/3 pP , ^ p2 Bp 



2(1-6) 



(23) 



(24) 



where g^^ is the metric tensor in the transverse components. We can further express P^j^ in terms of Mandelstam 
variables. 



stu 



(s + g2)2- 

Finally our result for the real corrections is given by 



(25) 



d(P^iAa(5x)) 



dxBdydzh 



Real 



as ( 4:TT^ 



2\ e 



V Q2 ; r(i-e)i 



-z-^l- z)-''x'il-xy 



x-^Tg,Fix,x)Dl+Tg^Fix,x)N^+Tq^Fix,XB)N^ 



(26) 
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where we have used the expression for the two-body phase space integral in n = 4 — 2e dimension 



12 / r(l-e) 
The hard-part coefficient functions in Eq. (PBl) are 
f s 



dzh I f ^■nV 1 



D 



« 2NcS + Q2 
1 f 



^ t s st 



« 2A^, M(s + Q2) 



[{S + Q2)3 + q2(^2 _ ^^2) ^ ^^^^ ^ _ ^(g2 ^ ^) ^(g2 ^ ^)(g2 ^ ^) ^ 2Q2s)] 



Cf + 

t + u 2Nc 



si(s + Q2) 



- Q^ii^ + diQ^ + u){i+u)] 



(27) 

(28) 
(29) 

(30) 



C. Combining the real and virtual corrections 



To combine real and virtual corrections and demonstrate how the divergences cancel out, we need to perform the 
e-expansion for the real corrections in Eq. (|26p . To proceed, we realize that 



Let us define the following common factor: 



1 - e 



X X 



z-'{l- z)~'x'{l-xy 



y\^e carry out the e-expansion for the / x {D^ ,N^,N'}) terms. Using the following formulas [2 



-e-l 



e ^ ' {l-z)^ 

— 5(1 - x) + - — 3_ 

z"'(l - zy = 1 - elnz - eln(l - z), 
a;'^(l — i)^' = 1 + eln.x — eln(l — x), 



r'-{i-z) 



x^{l-x)- 



ln(l - z 
I- z 
f ln(l — X 



\ l-x J ^ 1- 



In z 

'Y~z' 

Inx 



(31) 



(32) 



(33) 

(34) 

(35) 
(36) 



we present the expanded results for these hard-part coefficient functions. For the so-called derivative term Z?^ we have 
1 ' 



--{I + x^)5{l -z) + {l-z)+ ,f^X'^'''' - (1 + In 7^-5(1 - z) - 2x5{\ - z) 
t (1 — z)+ 1 — X 



(37) 



For the first divergent piece we further perform integration by parts to convert x-^Tq^p{x, x) to the function Tq^p{x, x) 
itself: 



1 f dx d ^ , ^ , 9, 1 f dx ^ , ^^9 
— - x-T,A^, x)il + .^) ^ ^ -TqA^, X) [2x^ - 2.(1 

Thus, the term associated with leads to the following expression: 
/ ? ("0 ^)^'''^^^'^)2^ ^^^^ ^ 



(38) 



+x-T,M^,x) — 



(1-5) + 



(l-x)^ +2xz 
(1-5) + 



(1 + In -5(1 -z)- 2x6(1 - z) 



1-x 



(39) 
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Likewise, we have 



1 



1 l + z" 



2iV, 



e(l-z)+ 



5(1 -x) 



1 2x^ -3x^-1 



{l-xU 



5{l-z)-25{l-x)5{l-z) 



+5{l-z) 



-(1 -x){l + 2x) In 



1 -X 



f \n{l-x) 
V 1 - X 



2(1 -i) + 2' 



Ini 
l-x 



(1 + z) lnz(l-z)--2 



2x-^ - - 1 



1 + z 



I xN: = 



(l-i)+(l-z)+ (1-^)+ 
1 



/ln(l - z) 
2(1 -i) I 



2z 



(1-z), 



In z 
I 

'l - z 



(40) 



+ 2iV, 
1 l + x 



4-5(1 - £),5(1 - z) + -5{l ~ x)5{l -z)~ - /^f" 8{l - x) 



e(l-z) + 



-7^-r^^(l - 5) + 25(1 - - 5) + ^^^.^ 

e(l-.x)+ (l-a;)+(l-z)+ 



+5(1 - z) 
+5(1 -i) 



In- 



\-x 



ln(l - x) 
l-x 



l + x 



1 — X (1 — x)^ 



-(1 + z) lnz(l - z) + 2 



/ln(l - z) 

V 1-5 



In; 



2z 



1-5 (1-5)h 



(41) 



Now let us combine the resuhs for the real corrections in Eqs. ([39]), (j40|) . and (j4T|) with the virtual correction in 
Eq. (fT5|) . First, for the double-pole 1/e^ term, which represents a soft-collinear divergence, we find that they cancel 
out between real and virtual diagrams. Specifically, the term ^ l/2Nc in the soft-pole contribution cancel the 
corresponding term in the hard-pole contribution A'^^, which leaves the remaining 1/e^ term in A^^' with a color factor 
Cf- This remaining term exactly cancels the 1/e^ term in the virtual diagrams. 

Now we turn our attention to the 1/e term. By adding the corresponding terms in real and virtual diagrams, we 
end up with the following expression: 



(42) 



5 -- [<5(1 - x)Tq^F{x, x)Pqq{z) + 5(1 - z)Pqg^qg ® Tq^F{x, Xx)] , 



where Pqq{z) is the usual quark-to-quark splitting kernel, as in Eq. (jH). Thus, the first term, the part containing i, 
is just the collinear QCD correction to the bare leading order fragmentation function I?^'^^(z/i) (after including the 
pre-factor (47r/xVQ^)7r[l - e]): 



Dq^u{zu) - D^q\^{zu) + (-|) Dq^,{z)Pqqiz), 



where we have adopted MS-scheme, and 

— — 7b + ln47r. 

i e 

On the other hand, the second term —\Pqg^qg ® Tq,F{x, xx) is given by 
1 



(43) 



(44) 



-- ) Pqg^qg ®Tq^F{x,Xx) 



1 

e 

Nc 
' 2 



^q^F(x, x)Cf 

l + x 



1 + i^ 3,,_, 

^ h -5(1 - x) 

{l-x)+ 2 ^ [ 



Nc6{l - x)Tq,F{x, x) 



[l-x). 



Tq^F{x,Xx) 



l + x' 
{l-xl 



-Tq,F{x,x) 



(45) 
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from which we immediately obtain the colhnear QCD correction to the leading order bare Qiu-Sterman function 



T'fp{xB,XB) as follows: 



l + x 



27r 



dx 

X 



Tq^F{x,x)Cf 



— + -8{l-x) 



Nc5{l - x)Tq^F{x, x) 



' 2 



-Tq^pix, XX) 



1 + i 



C.2 



-Tq^F{x,x) 



From this equation, we can obtain the evolution equation for the twist-3 Qiu-Sterman function as follows 



(46) 



51n^2 9 





^ dx^ 
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- X 






xY 



— + -5{l-x) 

(l-x)+ 2 



- Nc5{l ~ x)Tq^F{x,x,n ) 



-Tq^Fix,XX, ) - 



(l-x) 



-Tq^F{x,X,^I ) 



(47) 



22-24|. In particular, we are able to reproduce 
17H19|. This piece arises as follows: there is a 



This result agrees with earlier findings from different approaches [20i 
the term —NcTq^F{x,x, /j?), which was missing in some early works 
boundary term oc f 2^^(-'^ ~ from D^, as in Eqs. and and this term cancels the same term with opposite 
sign in Nq in Eq. PO)) . On the other hand, the hard-pole contribution Nq contains the following term: 



zCf 



2Nr 



-6(1 - x)S{l -z)^[Cf + ^] -S{1 - x)5a - z) 



-Nc5{l - x)5{l - z)] , (48) 



which gives exactly the ~NcTq,F{x,x, ^'^) term to the evolution of the Qiu-Sterman function Tq^F{x,x, fi^). 

Now, let us turn to the finite NLO corrections to the hard-part coefficient function. After MS subtraction of the 
collinear divergences into the fragmentation function Dq^h{z, IJ^'^) and the twist-3 Qiu-Sterman function Tq^Fix, x, /i^), 
we obtain the NLO correction for both soft-pole and hard-pole contributions to the P/jj^-weighted transverse spin- 
dependent cross section: 

"1 (j^x dz 

— Tq^F{x, X, ^l'^)Dq^h{z, fl'^)S{l - X)S{1 - z) 



d{Phj_Aa{S±)) ZhCTo 



dx B dydzh 



2 2 



+5{l - z)Pqg^qg (g) Tq^F{x, XX, ^^)] 

1-z {l-x)'^ + 2xz 



[(5(1 - x)Tq^F{x,X,fl'^)Pqq{z) 



+ X^TqMx^X,^,')^ 



+Tq^F{x,x, ^'^)6{l - z) 



1 



z{l-z)+ 
{2x^ - £ - 1) In 



(5(1 - z) (1 +i:^)ln 



2x 



1 - X 



- 2 



/ln(l -; 



l-x 
2x{2 ~ x) 



In; 



+Tq^F{x, X, Ai^)^(l - x)Cf 



-(l + z)lnz(l-z) + 2 



2z 



(1-^)h 



+Tq^F{x,X,^l'^) 



1 



2NcZ 



2x^ - 3^2 - 1 



l + z 



{l-x) + {l-z)+ (1-.t)_ 



+Tq^F{x, XX, /i^)(5(l — z) 



+Tq^F{x,XX, f/) 



2 

1 -1- xz^ 



In- 



{l-x)+{l-z)+ 



l-x 
Cf 



ln(l - x) 



l-x 



V l-x 

ln(l - z) 
1-z 

2(1 -X) 

2 In a; l + x 
_^ l-x (1 - i:)+ 



l-x 



In z 
I 

'l-z 



1 



2NcZ 



Tq,F{x, X, /i^)6CF(5(l — x)5{l — z) 



(49) 



Just like the NLO correction to the spin- averaged cross section in Eq. ([7]), the logarithms containing the factorization 
scale together with the splitting functions determine the evolution of the twist-3 Qiu-Sterman function and the usual 
unpolarized quark-to-hadron fragmentation function. Eq. (149^ is the main result of our paper: once combined with 
the NLO spin-averaged cross section in Eq. ([7]), one will be able to compute the P/i^-weighted Sivers asymmetry as 
defined in Eq. 
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IV. CONCLUSIONS 



We calculated the next-to-leading order perturbative QCD corrections to the transverse momentum-weighted Sivers 
asymmetry in semi-inclusive hadron production in lepton-proton deep inelastic scattering. Specifically, we demon- 
strated in detail how to evaluate at NLO the P/i_L-weighted transverse spin-dependent differential cross section. We 
found that the result can be written as a convolution of a twist-3 quark-gluon correlation function (often referred as 
Qiu-Sterman function), the usual unpolarized fragmentation function and the hard coefficient function. In the course 
of this calculation we showed that the soft divergences cancel out between real and virtual contributions, and that 
the collinear divergences can be absorbed into the NLO twist-3 Qiu-Sterman function of the transversely polarized 
proton and the unpolarized quark-to-hadron fragmentation function. Such a procedure also provides an alternative 
way to identify the evolution equation for the twist-3 Qiu-Sterman function. We found that our evolution equation 
agrees with those derived previously from different approaches. Using our NLO results for both the spin-averaged and 
-P/i±-weighted transverse spin-dependent differential cross section, we plan to study in the future the P/ij^-weighted 
Sivers asymmetry. We anticipate that our findings will have important phenomenological applications relevant to the 
experimental programs at HERMES, COMPASS, and Jefferson Lab. 
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